% % % cell_num=[3 7 12 13 15 19 22 23 24 26 31 43];
% % % activity_matrix={};
% % % cnt=0;
% % % for j=1:length(cell_num)
% % %     cnt=cnt+1;
% % %     act=ftfil(ROIs(cell_num(j)).activity,18,0,10);
% % %     act=psmooth(act);
% % %     raw_act=act/median(act);
% % %     activity_matrix{cnt}=raw_act;
% % % end

activity_matrix2={};
for j=1:length(cell_num)
    activity_matrix2{j}=activity_matrix(j,:);
end
frame_times=sess.frame_times;
aux_data=sess.aux_data;
stim_id=round(aux_data(12,frame_times)*2);
run_play_id=round(aux_data(2,frame_times));
run_id=round(aux_data(9,frame_times));

% % % figure;plot(run_id)

max_act=[];
min_act=[];
for g=1:length(activity_matrix2)
    max_act(g)=max(activity_matrix2{g});
    min_act(g)=min(activity_matrix2{g});
end

maxact=ceil(max(max_act));
minact=floor(min(min_act));

% % % obj_pos=(9/18);
% % % max_pos=max(run_id);
% % % min_pos=min(run_id);
% % % maxP=max_pos+10;
% % % minP=min_pos+10;
% % % rang=maxP-minP;
% % % pos=rang*obj_pos;
% % % obj_pos2=(minP+pos)-10;
% % % 
% % % el_pos=[];
% % % dire=0;
% % % stop=1;
% % % for u=1:length(dir_ind_ft)
% % %     if u>1
% % %         if dir_ind_ft(u-1)~=dir_ind_ft(u)
% % %             stop=1;
% % %         end
% % %     end
% % %     if dire==0
% % %         if u==1
% % %             if run_id(u)<obj_pos2
% % %                 go=1;
% % %             else
% % %                 go=0;
% % %             end
% % %             if go==0
% % %             elseif run_id(u)>obj_pos2
% % %                 if stop==1
% % %                     el_pos(u)=u;
% % %                     stop=0;
% % %                     dire=1;
% % %                 end
% % %             end
% % %         else
% % %             if run_id(u)>obj_pos2
% % %                 if stop==1
% % %                     el_pos(u)=u;
% % %                     stop=0;
% % %                     dire=1;
% % %                 end
% % %             end
% % %         end
% % %     else
% % %         el_pos(u)=0;
% % %     end
% % %     if dire==1
% % %         if u==1
% % %             if run_id(u)>obj_pos2
% % %                 go=1;
% % %             else
% % %                 go=0;
% % %             end
% % %             if go==0
% % %             elseif run_id(u)<obj_pos2
% % %                 if stop==1
% % %                     el_pos(u)=u;
% % %                     stop=0;
% % %                     dire=0;
% % %                 end
% % %             end
% % %         else
% % %             if run_id(u)<obj_pos2
% % %                 if stop==1
% % %                     el_pos(u)=u;
% % %                     stop=0;
% % %                     dire=0;
% % %                 end
% % %             end
% % %         end
% % %     else
% % %         el_pos(u)=0;
% % %     end
% % % end

size=(length(activity_matrix2)+3)*7.5;
width=length(frame_times);
ids_a=unique(stim_id);
ids=ids_a(abs(unique(stim_id))>0);

figure('color','black')

axes('position',[0.01 0 0.98 1])
hold on
%title(regexprep(mouse_id,'_',' '),'color','r')
axis([-500 width -1 size])
axis off

for h=1:length(ids)
    area(1:length(stim_id),size*(stim_id==ids(h)),'facecolor',[0.2 0.2 0.2]*(1-(h/15)),'linestyle','none')
% % %     text(-500,ids(h),num2str(ids(h)),'color','yellow')
% % %     text(-200,ids(hnd),'-','FontWeight','bold','color','blue')
end

% % % for t=1:length(el_pos)
% % %     if ~el_pos(t)
% % %     else
% % %         area([el_pos(t) el_pos(t)+79],[size size],'facecolor',[0.2 0.2 0.2],'linestyle','none')
% % %     end
% % % end

plot(stim_id,'color','yellow')
% % % plot(run_id,'color','yellow')
for n=1:length(activity_matrix2)
    hold on
    plot(activity_matrix2{n}+(size-(n*maxact)),'color','red')
    text(-500,(size-(n*maxact)),num2str(cell_num(n)),'FontWeight','bold','color','red')
% % %     if n==length(activity_matrix2)
% % %         
% % %         axes('position',[0.01 0.01 0.98 0.04])
% % %         hold on
% % %         axis([-500 width -5 5])
% % %         plot(run_id,'color','white')
% % %         axis off
% % %         
% % %     end
end
axes('position',[0.01 0.05 0.98 0.04])
hold on
axis([-500 width -5 5])
plot(run_play_id,'color','blue')
plot(run_id,'color','white')
axis off

set(gcf, 'InvertHardCopy', 'off');
filename=strcat('C:\Users\zmarz\Desktop\LEC_results\',curr.mouse_id,'\',curr.date_id,'\','plot');
print(gcf,'-dbmp',filename)